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We compare fully differential exact results for the virtual photon correction to single hard photon 
• bremsstrahlung obtained using independent calculations, both for e + e~ annihilation at high-energy 

colliders and for radiative return applications. The results are compared using Monte Carlo eval- 
uations of the matrix elements as well as by direct analytical evaluation of certain critical limits. 
Special attention is given to the issues of numerical stability and the treatment of finite-mass cor- 
J_h ' rections. It is found that agreement on the order of 10 -5 or better is obtained over most of the 

range of hard photon energies, at CMS energies relevant to both high energy collisions and radiative 
return experiments. 
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The single hard bremsstrahlung process is critical in precision studies of Standard Model physics 
over a wide range of scales, from the 1 GeV regime, where it is exploited in radiative return studies 
of intermediate-energy hadron physics, to the 100 GeV regime and beyond, where it is used for 
' radiative return studies of the Z resonance. Construction of a future linear collider would push the 

energy range of interest toward the TeV regime. Precision studies of per mille level effects in any of 
these regimes requires exact calculations of all QED effects at order a 2 . For radiative return, this 
means that the virtual photon corrections to single hard photon radiation must be included. 

43 ■ 

Exact results for the virtual corrections to single hard bremsstrahlung have been studied by 
a number of groups. P, 0, 0, 0] Comparisons of these results are important to be certain of their 
correctness at a high level of precision. In ref. 3|, we have previously presented comparisons of the 
results in refs. [HQ. El- an d m general a very good agreement was found. However, it is clear from 
the comparisons in ref. that small differences appear at the level of the NNLL (next-to-next-to- 
leading- log) , suggesting differences in the levels of "exactness" in the calculations. 

These differences are largely due to how differential the results are, and how the mass corrections 
are included. Both ref. !3| (JMWY) and ref. pj include mass corrections, but the former is fully 
differential, while the latter averages over the photon angle. The results of ref. are fully differential, 
but the mass corrections are neglected. These comparisons therefore are not a complete test of the 
results in ref. Q at the NNLL level. 

Another calculation of the exact virtual corrections to single hard bremsstrahlung has appeared 
recently. An exact differential result including mass corrections needed for an accurate calculation 
of collinear emission is presented in ref. Q (KR). This result was obtained for use in radiative return 
studies of hadronic final states. 0, ||| QJ The virtual correction without explicit mass termsQ was 
obtained for use in tagged photon experiments, and this was extended to the small angle regime for 
untagged experiments with the inclusion of finite mass corrections Q, which become significant in 
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collinear emission, even in very high-energy collisions. 

The results of KR have been implemented in the PHOKHARA MC [1 B S E3, while the 
results of ref. Q have been implemented in the JCJC MC ■ In re f- > an indirect comparison of 
the JMWY and KR results has been reported via comparison of the PHOKHARA and K.IC MC's. 
Agreement at the per mille level was reported for selected observables. Further comparisons of 
these MC's were reported in ref. where the effect of including YFS (Yennie-Frautschi-Suura) 
exponentiation 14] in the K.JC MC versus the unexponentiated approach of PHOKHARA plays an 
important role in the differences found. 

Presently, our interest is not in comparing the MC's but in comparing the order a 2 matrix 
elements for hard photon emission used in each. A direct comparison of the matrix elements JMWY 
and KR was reported in ref. [l]|, using a common MC evaluation of both expressions. It was found 
that the expressions agree to 1.5 x 10~ 5 units of the Born cross section for a photon energy below 
95% of the beam energy at 200 GeV CMS energy. 

This level of agreement was the result of careful analytical work to cancel large terms in the 
expressions of KR. The published matrix elements of KR and JMWY have very different forms, 
which is evident in part from the appearance of mass corrections proportional to (pi ■ k)~ 2 and 
{pi ■ fc)~ 3 in the expressions of KR, where pi is an incoming fermion momentum and k is the emitted 
hard photon momentum. Such terms are absent in the JMWY expression, and we have verified that 
they cancel exactly in the KR expression as well, although this is not immediately obvious. This 
cancellation should be implemented analytically to obtain a stable MC comparison of the results. 
Such an analytical cancellation is presented here. 

The functions used to construct the stable comparisons were introduced in ref. 0| . A proof that 
the massless limits of the KR and JMWY expressions agree at NLL (next-to-leading-log ) order was 
sketched and the size of the difference in the mass corrections was calculated in ref. ■ MC results 
for a comparison of the matrix elements in the 1 GeV radiative return regime were also obtained. |18| 
In this paper, we will present details of the stabilization procedure needed for very high precision 
comparisons, and present new MC comparisons which use these versions of the matrix elements, 
which have been improved since the publication of refs. [IE EE E3, E3 

for better behavior at lower 

CMS energies. 

An important part of this comparison is an examination of the finite-mass corrections used in the 
virtual corrections of KR and JMWY. The electron mass terms containing denominators of the form 
Pi ■ k give contributions of the same order as some of the "massless" contributions upon integration. 
Therefore, the finite mass corrections of both calculations are compared explicitly, and it is shown 
that both mass corrections are identical to next-to-leading-log (NLL) order. This supplements the 
result that the massless results agree to NLL order, showing that the complete virtual corrections 
agree up to NNLL terms. These NNLL terms are the subject of the numerical comparisons, both 
with and without mass terms. 

We need to point out that both approaches to the exact mass corrections discussed here, that 
used by JMWY and that used by KR, should have the property that, upon integrating over the 
final phase space, the resulting integrated cross section should contain all finite terms that do not 
vanish for m 2 /s — > 0. Both approaches use expansions therefore in the parameter m 2 /s. The 
JMWY approach uses the expansion constructed by ref. 0] . The KR approach uses the expansion 
of the respective leptonic tensor coefficients through terms of order m\sj(2k ■ Pi) 3 , where k is the 
radiated real photon four momentum and pi is an incoming e~ or e + four momentum. It is therefore 
important to know the relationship between these expansions both numerically and theoretically. 
We address this relationship matter in what follows. 

The structure of this paper is as follows. In section [H] the matrix elements of JMWY and 
KR are displayed explicitly in a compatible format appropriate for comparisons, without the mass 
corrections in each case. Both forms of the matrix elements are shown in a form explicitly stable in 
collinear limits, which is suitable both for MC evaluation at very high precision, and for an analytical 
investigation of the collinear limits. A comparison of the collinear limits is presented in section ITTT1 
demonstrating that both matrix elements agree to NLL order when the explicit mass terms are 
omitted. Section llVI examines the mass terms which are added in each case, and verifies that these 
both agree as well to NLL order, and to NNLL order in the soft limit. A numerical comparison of 
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the matrix elements is presented in section both at 1 GeV and 500 GeV CMS energies. Finally, 
section IVTl summarizes our conclusions, and the Appendix presents details on the special functions 
introduced to stabilize the collinear limits. 



II. VIRTUAL CORRECTIONS TO INITIAL STATE RADIATION 



The process we consider is initial-state radiation (ISR) in e + e~ — > fi + fi~ , where one real and 
one virtual photon are radiated on the electron-positron line. In ref. Ej, the matrix element was 
evaluated using helicity amplitude techniques of CALKUL[2(J, Xu et al 21] and Kleiss-Stirling p^ . 
together with algebraic simplification using FORM [23|. The mass corrections were then added 
following the methods in ref. Hi, after checking that the exact expression for the mass corrections 
differs from the result obtained by the latter methods by terms which vanish as m 2 /s — > in the 
integrated cross section, where m e is the electron mass and s is the squared CMS energy. 

In ref. [j], the same ISR matrix element's interference with the Born process is calculated in 
terms of a leptonic tensor in such a way that all mass effects relevant in the collinear limits are 
obtained using an expansion in powers of n%/pi ■ k. Thus, comparison of the two sets of results 
gives important information on the two different methods of calculation and on the two different 
treatments of the mass corrections. 

The JMWY result used in these comparisons will be evaluated without Z boson exchange, to 
match the calculation of KR. We denote the four momenta and helicity of the e~ , e + , /, and / as 
pj and Xj, j = 1, ...,4, respectively. The tree- level cross section for single hard ISR for QED alone 
may be written as 



da 



ISR(O) 



dQdri dr 2 



16(4tt) 4 E 



(1) 



where the squared matrix element, summed over helicities, is 



ISR(O) 



1(!e6 r (t?+o(i 2mln 



s 2 s'rir 2 



s'r 2 



2m 2 e r 2 
s'ri 



(2) 



with s = (pi+p 2 ) 2 , s' = (P3+P4) 2 , h = (P1-P3) 2 , t 2 = O2-P4) 2 , «i = {pi-pa) 2 , u 2 = {p 2 -pz) 2 , 
ri = 2pi ■ k/s, and the solid angle d£l refers to the direction of p% — pa relative to p\ — p 2 in the CM 
frame of the outgoing fermions. The explicit mass corrections in eq. (J2J are those obtained using 
the method of ref. 19] . 

The same tree-level matrix element may be expressed in terms of the leptonic tensor L^ v of 
refs. 0, 0- I n the same notation, the leading-order contribution to the leptonic tensor may be 
written (in the normalization of ref. 0) 



L - 



s(s') 2 rir 2 



\ a 00 sr l i a ll Pi Pi + a 22 P2P2 ' a 12 VP1P2 + P2P1 



n} 



(3) 



with 



where z — s'/s 
muon line, 



(0) 

'on 



(l-r 1 ) 2 -(l 
-4z+ 



2m 2 zv 2 



sr\r 2 



= -Az- 



8m 2 e r 2 (o) 



-8mi 



L2 



(4) 



sr 2 sri 
v with v = ?*! + r 2 . This may be contracted with a final state tensor for a 



M"" = e 2 {p&l + pfo% ~ (P3 • P4)tT } , 



(•5) 
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to obtain the same squared matrix element as in eq. J3J, 
£ |.Mi SR(0) 2 = WL$M» V 

\ -Ss'(2(Ioq + ) + ttfi tlUl + a 22 ^ 2 U 2 + a 12 (*1* 2 + U1U2) } , 

s(s')^rir 2 L J 

(6) 

where the coefhcient is chosen to match the normalization of eq. (J2J), and explicit mass terms are 
kept only when enhanced by collinear factors. The mass terms in eq. J2Jl and eq. 10 agree in the 
collinear limits where is of order to^/s. 

The initial-state virtual photon corrections to the cross section may be expressed as 

where the matrix element for hard photon initial-state emission with one virtual photon may be 
expressed as 

Mf R U = ^Mf R V\f + Ah + hh), (8) 
where fi are scalar form factors and li are spinor factors defined in ref . Q : 



h = a\ 3 8x 1 {pi,k)s^x 1 (p 2l k) 

S\ 3 (p4,P2)s-\ 3 (p2,P3,) ~ SA 3 (P4,Pl)S-A 3 (pi,P3) 



S_ CT (p 1 ,p 2 )S- (T (p3,P4)sJ(F21,P34) 
S Xl (pi, fc)s_ Al (P2j fc)sA 3 (P4, fc)s^A 3 (P3 , k) 
S-a(Pl,P2)S-a(P3,P4)sl(Pl2,P34,) 



(9a) 



7 2 — A1A3 r r , (9DJ 



where is the helicity of the fermion with momentum pi and a is the photon helicity. The spinor 
product is s\(p,q) = u_\(p)u\(q), and 

The form factors fi are given (without mass corrections) by 

f = 4nB YFS ( S ) + 2 (L - 1 - ztt) + + r ^+Il> j m + Z7r \ 

1 — r 2 (1 — ri)(l — r 2 ) L \zJ > 

3w + } Lfi (-«) + - , i2i(r 1( r 2 ) + r 2j Ri(r 2) n), (11a) 

1 - r 2 J (1 - r 2 ) 

A = 2(1 r 2 ) + (1 - rxKl - r 2 ) + 2^)) fa (?) + 

1 „ , s „ , si V 



, --Ri(ri,r 2 )+r 2 ii 2 (ri,r 2 ) ^ + - {fli(n, r 2 )<5 ff ,i + Ri(r 2 , r^S^-t} , (lib) 
1 - r 2 2 J 4 

/2 = V+ r- - 2vzLf 3 (-V 

(l-ri)(l-r 2 ) 



z 



(l-ri)(l-r 2 ) V 2(1 -ri)/ I V z 

z f 1 



l-r 2 12 
ri - r 2 



Hi(ri,r 2 ) + (2-r 2 )ii 2 (ri,r 2 ) 
{i^i(ri,r 2 )<5 (7> i+ J R 1 (r 2 ,ri)5 (7 _i} (11c) 
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for photon helicity a = X%, When a — — Ai, ri and r2 must be interchanged. Here, L = hi(s/m 2 ) 
is the "large logarithm" which may be used in a leading log expansion of the results. Thus, leading 
log (LL) refers to order L 2 , NLL to order L, and NNLL to order L° = 1 in this expansion. 

The infrared divergence is contained in the virtual YFS form factor defined by [l4| 



2 



(mz \ „ An' 

2 In — 1 + 1 (L - 1 - m) - L 2 - 1 + — + in (L - 1) . (12) 
mi / 3 



We also make use of functions 



Ri(x,y) = Uti-x) { In ( — — ) ~ 2m 



2(1 - X ~ V) Sfx f X(1 '"' yr ) , (13a) 



1 — X \ 1 — X 1 — X 

J? 2 (x,y) = U 1 (-x)-hi 2 1 (-x) + ^-U 2 (-x)+ f -Lf 2 

2 x + y [x + y){l-x-y) 

+ Lf.^dn. + ^+fi-^) 2 Sfjg ^ 1 -«->) ), (13b) 
\ 1 — x y \l — x 1 — x / 

with Lf n (x), Sf n (a;,?;) defined in eq. (|Ala|) of the Appendix. 

The functions Lf„ and Sf„ play an important role in canceling large factors in differences between 
logarithms or dilogarithms with slightly different arguments. Therefore, although the expressions 
for /o,/i,/a are analytically equivalent to those of ref. Q, they are more stable when evaluated 
numerically. The functions R\ and i?2 are related to the function R used in ref. 3] by 

Ri(x,y) = -R(x,y), 
x 

M*,V) = -llRi(x,y)-U 1 (-x)+ln(-^ (14) 



The functions Ri are defined so that both are finite for x,y — > with x,y positive. This implies 
that fx and f% are finite for r%,T2 — > 0. The spinor products I\ and I2 vanish in this limit, since 
\s\(j>i, k)\ 2 = sri. Therefore, the f\ and f 2 terms are absent in the collinear limits. 

The virtual corrections aff to the coefficients a,j are calculated without explicit mass terms 
in ref. 8J. The infrared divergence in that result is canceled by adding the contribution of an 
additional soft real photon with energy below a cutoff u m i n as a fraction of the beam energy. To 
facilitate comparisons, the additional soft photon contribution will be removed, so that the pure 
virtual corrections are compared in each case. This requires decomposing the infrared-divergent 
term 



a\f = ag 3 \2{L - 1) In v mia + I (L + In z) - 2 + ^ \ 1 1 :, ) 



in ref. 8] into a virtual contribution afj -Fir and a contribution due to an additional real soft photon. 
The correct decomposition is found to be 



a]f = a$(F m + 2irB YFS ( s ,v mhl )) (16) 



with 

F m = 27rRe B Y fs {s)+L-l + ~biz (17) 
and the real YFS soft photon form factor 

f to 2 1 1 7T 2 

2nB YFS { Sl v min ) = (£- 1) |ln-f + 21nw min | + -L 2 - —. (18) 
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This can be verified by comparing the soft limit of the virtual correction to the known result, as is 
done at the end of this section. 

We will use the functions Lf „ and Sf „ defined in the Appendix to stabilize the results and clarify 
the collinear behavior. Then ay = a\j + 7a with 

Ji) _ 01 f„(P) „ (1,0) „(i,m)\ h q\ 

where the massless parts a,^' ^ of the non-IR-divergent part of the virtual corrections are given by 

+ {n + ?[ 1 + (i-^]}{,.n©L tl «-n)-S fl (-Il.- ? )} 

+ (ri^r 2 ), (20a) 
a{\ 0) = 3 + z+ -61nz-4Lfi(-u) 



- 4(1 + n - nr 2 ) {Lf x (-t>) + Lf 2 (-v)} + 2 Q - z(l - ra)) Lf x ( 
+ 2(1 +*)((!- r 2 ) Lfi(-ra) - Lf x f— ^) - - Lf t f ^ 



nr 2 



1 — r 2 / z V z 
+ 2 Lfi(-n) - I Lf! (^) + 2r 2 zln (^-) {Lfi(-r 2 ) + Lf 3 (-r 2 )} 

+ ri«ln(^) in (jz^) +4Lf 1 (-r 1 ) + 2(l-r 2 ) 2 Lf 3 (-r 1 )| 



2r 2 Sfi I — ^, + - Sf 2 



z z 1 z \ z z 
- 2(2 + n - 2ra)Sf! ( ^ + 2z Sf 2 ( , (20b) 

4 2 0) = ali 0) (ri-r 2 ), (20c) 

(1.0) z(l + z) V 

a 12 = — — 7 + 2r l r 2 



2(l-n)(l-r 2 ) 2 
- z Lfi(-w) - 2z(l - nra) Lf 2 (-«) - z Lf x (^) 

+ 2z(nr 2 + r 2 - 1) I y-i— + Lf 2 (-n)| In (y) 



ri 

Infra) _ 2 (l-r 1 )Sf 1 (^,^U2Sf/-' ! 1 



(1 — r 2 ) 2 \ z J \ z z ) \ z z 

+ (n «-> ra). (20d) 

These expressions are exactly equal to those of ref. Q, but written in terms of the functions Lf„ 
and Sf n . 

A comparison of the collinear limits of the two matrix elements will be presented in the following 
section. A numerical study comparing these two forms of the virtual photon correction to hard 
bremsstrahlung will be presented in section [V] together with a comparison of the mass corrections 
added to obtain the complete matrix elements in each case. 
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III. NLL COMPARISION OF THE MATRIX ELEMENTS 



In this section, we will verify that the expressions of ref. Q and ref. agree analytically to 
NLL order, meaning that they agree through order L when integrated over n. For this purpose, 
we may examine both expressions in the limit of collinear emission, which is sufficient to determine 
the NLL result. We will consider the particular limit where r\ becomes small, so that the photon is 
collinear with the incoming electron. 

In the collinear limits, the virtual correction to the squared matrix element is proportional to 
the tree-level result, 

£ 2 Re (M\ SR(0) *M\ SR(1) ) = £ Re (/ NLL ) £ |^i SR(0) | 2 , (21) 

where (/(^ LL ) is the average of /□ over helicities, evaluated for small r\ or r 2 small, 

Re (/ NLL > = 27rRe B YFS {s) + L - 1 + 3 ln(l - n) + 2 Inn ln(l - r 2 ) 

- ln 2 (l - ri ) + 2S P (n) + ^ { }~ ri \ 2 + (n -> r 2 ). (22) 

1 + (1 - ny 

Since the /i, / 2 terms in eq. (jSJ) both vanish in collinear limits, they do not contribute to NLL 
order. In the limit when r% — > 0, the NLL expression may be written in the form 

Re (/ NLL (n ->())} = F m + 2 In z In n - In 2 z + 2 Sp (v) + — ^ . (23) 

_1_ — \~ z 



The result of ref. 8] can be compared to this expression using eq. JBJ with af^ replaced by 
in the collinear limit. The virtual correction to the squared matrix element may be written 



£2Rc (Mf^'MfW) =16L^M^ = ^2F m y £\ M f m 



+ 



s(s') 2 rir 2 



2a 



(1,0) 



(1,0), 

ail 



,(1,0) 



t 2 U 2 + Oi2°'(*l*2 + U\U2) 



(24) 



In the collinear limit, this expression should again be proportional to the tree level result, eq. @ 
without the mass terms. When r\ is small, ti = zt\ +0(r\) and u 2 — zu\ +0(r\), allowing eq. (|24() 
to be written as 

£ 2 Re (^r ) )^{/'Ehr , f w 

Aj,<7 I Aj,<7 J 

with 

a oo v 4/ 
A = ^+za^-2a^-Aa&°\ (26) 

where clqq = — (1 + z 2 ) is the massless limit of the tree- level coefficient. When r 2 is small, the same 
form is obtained with an and a 22 interchanged. 
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We may check that A = and /' = Re fo in the collinear limits. Evaluating the coefficient 
functions in the limit when k is collinear with pi gives 

4o 0) (n^0) = -y + (f + z 2 )|iln 2 z-lnzlnn-Sp(«)|, (27a) 

4i'° 3 (n ->0) = 3 + 2z- 21nn |zlnz+ U^-v) + U 2 (-v) - ^(1 + v) 
+ In z {Lf i (- v) + 2 Lf 2 (- u) - 8} + z In 2 z - 2zSp (y) 

- 11 Lfi(-u)-Lf 2 (-w)-5Lf 2 (-w) + 2Sfi(0,w)-2Sf 2 (0,i;), (27b) 
a^ 0) (n ->0) = i + 2z + lnr 1 |i-4zlnz-2Lf 1 (-v)-2Lf 2 (-u) 

3 .1 . „ , 2 , /'2 



+ lnz|Lfi(-w) + 2Lf 2 (-w)---4|+2zln z-Lfi(-u) I - + 7 

- 4zSp(w)- Lf 2 (-u)-5Lf 2 (--y) + 2Sfi(0,'y)-2Sf 2 (0,f), (27c) 

a i2° ) ( r i- >0 ) = 3 + z + ]nri{l-2Lfi(-u)-2Lf 2 (-«)} 

+ hxz{Hi(-v) + 2 Lf 2 (-«) - 6 - 22} - 10 Lfi(-w) 

- Lf 2 (-w)-5Lf 2 (-v) + 2Sfi(0,-y)-2Sf 2 (0,i;). (27d) 

Combining these relations confirms that A = and /' = Re fo, so the results have the same 
collinear limits without explicit mass terms. This implies that the results without mass terms agree 
to NLL order. 

Since the collinear limits agree, the soft limits also agree. Specifically, for v — ► 0, 

a 00 - U ' - fl 22 - a 12 - l 28 ) 

and /' = Re fo — 4-7T Re£?YFs( s ) + 2(L — 1). This verifies that the virtual part of the infrared 
contribution was identified correctly in eq. 



IV. COMPARISON OF MASS CORRECTIONS 



Explicit mass corrections are added in the JMWY result via the prescription of ref. and 
sec. 3 of ref. |2^|, leading to a correction which may be expressed as 

£2 Re = -4° M ^\$+g + i+g\ i ( 29 ) 

where the non-radiative virtual correction factor is 

/nr(s') = 47rBYFs(a') + 2 In (-O - 2 (30) 



with 

4ttSyfs(s') = 47T-B Y Fs(s)+hi2 [ 2 In - 2L -Mtt ) - In 2 z. (31) 

The first term in brackets in eq. (|29|l is due to e + line emission, and the second is due to e~ line 
emission. 

This mass prescription is designed to produce the essential mass corrections in the collinear 
limits when the photon is emitted along the electron or positron line. In these limits, = m^v/s, 
where i labels the momentum of the collinear incoming fermion line. The mass corrections of ref. 
may be compared analytically in this limit. 
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The explicit mass correction of KR is given by 

£ 2 Re (M\ SR{0) *M\ SR{0) ) =16L£T>M' u ', (32) 

where only the mass corrections to the virtual part of the leptonic tensor are included, so that it is 
calculated using the coefficient functions 



(to) 



^■ m) F m . + 4 m) ) (33) 



with Fir as in eq. ijTTI) and a|" : " 1 ' given by the explicit mass terms in eq. The functions aj*' m ' 

are given in ref. and may be written in terms of pi = sri/m 2 e and the stabilized functions Lf„ 
and Sf„ as 

a g, m ) = zr2 r_ ln;2 ,(2£-4hipi+lnz)+4Sp(u)}-^r 2 Sfi(l,-pi) 
Pi 

- ^zr 2 Ni (pi, ^ - 4 J + vr 2 N 3 (p-i) + {rx,p x <-> r 2 ,P2), (34a) 

Jim) = l^r lnz (2L-41np 2 + lnz)+4Sp(i;)}-4rif- + 2 N ) Sfi(l,-pi) 

+ ^ ^zN 3 (p 2 ) + (2z - -^j Ni {^2, 2 ^~2*2 ) + 6N M + 2N i(P2 

+ ^{4N 3 (p 1 )-N 1 (p 1 ,0)+2z(N 2 (p 1 )+N 4 (p 1 ))}, (34b) 
v 

«22 = a^'" l) (ri,pi r 2 ,p 2 ), (34c) 



M , — {4JV 3 (pi) + 47V 2 ( Pl ) + 2JV 4 (pi) + 2 Sfi(l, - Pl )} - r 2 iV! ( Pl , -) 



+ <-> r 2 ,/9 2 ). (34d) 

where the functions ATj are defined by 



ATi(p,w) = Lf 2 (p-l)- W Lf 2 (p-l), (35a) 

i\T 2 (p) = Sh(z,v-p)+U 1 (-v)U 1 (^-l 

\v 

= Sh(v,p-v)--]npH 1 ( V —^) (p<l), (35b) 



z \ z 

N 3 (p) = i(l-Inp-Sfi(l,-p)) 
P 

= -lnpU 2 (-p)- Sf 2 (0,p) (p<l), (35c) 
7V 4 (p) = Sf 2 (z ) u-/>)--Lf 1 (-7;)Lf 2 (^-l)-— 

= -Sf 2 (v, p-v)-^]npU 2 (j^j (P<!)- (35d) 

All of the expressions Ni are finite in the collinear limit p — > u: 

JVi («,«;) = Lfi(-«) Q + w^ , (36a) 

JV 2 (») = Lf Lfx(-z), (36b) 

AT 3 («) = -Sf a (0,«)-ln«Lf a (-w), (36c) 

- -i|Lf 2 (-«) + i(l + Lf 1 (-z))|. (36d) 

Also, Sfi(l, — /?) — ► Sfi(0,?;) — lnu Lfi(— v) in the collinear limit. All of the functions N vanish in 
the limit p — > oo, where masses are completely neglected. 
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In the collinear limit p± — > v, corresponding to e + line emission, the total mass correction 

(rn 
ij 



functions a ™^ (including the infrared part) simplify to 



a>™> = — ^\2zF m +(3z-l)(lnvlnz + Sp(v))-2zLlnz Inv 

STi I Z 



z In 2 z 



A m ) 



~}, (37a) 
- jlni; Q + 2z + 4z Lf 2 (-w)) + 1 + z 



sri 1 \z 



- 2z 2 Lfi(-i;) + z 2 Lf 2 (-i;)+4zSf 2 (0,tO}, (37b) 

4 = — ^- \ 8^ir -SLlnz- -^(1 + z-Az 2 - 8z 3 ) Inv + 81nulnz + 8Sp (v) 
sr\ ^ z 6 

+ 4z\nv (Lfi(-w) -Lf 2 (-u)) -41n 2 z + 6 Lfi(-v) -Lf 2 (-v) 

- 4z(Sfi(0,i;) + Sf 2 (0,i;))-^(l + z-2z 2 )j, (37c) 

Am) 
'12 



ST l 



2z\nv(U 1 (-v) - 2 Lf 2 (-u)) - 4 + — )lnv 



+ 4zLf 1 (-w)-zLf 2 (-w)-2zSfi(0,w)-4zSf 2 (0,t;)-i|. (37d) 

Note that there are no terms proportional to m 4 /r 2 , so that all such factors in ref. Q cancel when 
the terms are combined into the stabilized functions used here. 

To verify that this result agrees with the JMWY result to NLL order, it is necessary to check 
all terms with explicit logarithms L in the collinear limit. In the e + line emission case, where k is 
collinear with pi , 

£ 2 Re (Mr'^r) = ^^\-(4 + ul)f' m -zs>A m \ (38) 



with 



(m) (rn) (rn) 

fm = a_^_ + z_a^_a = ^.^^ 



4z 



,2 



2L In— 1 + 1 -biz -L 2 + O(L ) (39) 



™9 



and 



Therefore, 



Am) 



+ za^-4a^=O(L ). (40) 



to NLL order. This agrees with the e + line emission part of eq. (|29fl to the same order, since 

o 

777- 

/nr(s') =2Lln^--i 2 + 2L(l-lnz) + O(L ). (42) 

Therefore, the both expressions for the explicit mass terms agree analytically through NLL order. 

We can also check that the mass terms agree exactly in the soft collinear limit. Taking v — » 
in the collinear results for the coefficients a ** 1 gives 

a^=2F m , eft? =-1/2, o&°=8Fir + 3/2, a { ™ } = 1/2, (43) 

which gives to eq. (|38|l with f' m = 2Fr and A m = 0. In the soft limit, /nr. = Fir, so the mass 
corrections are identical in the soft collinear limit. 
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FIG. 1: The virtual photon contribution to the /3j residual for muon pair production, integrated up to a cut 
»„,„. Case (a) is for a CMS energy of 1 GeV, and case (b) is for a CMS energy of 500 GeV. 



V. MONTE CARLO COMPARISONS 



Since both of the results compared are fully differential, they are well suited to MC implemen- 
tation. The JMWY results were developed for the JCJC MC described in ref. ^lj, which is for high- 
energy fermion pair production. The results from ref. 4] are implemented in the PHOKHARA MC 
developed for radiative return applications at $ and B factories. Structurally, these MC's are con- 
siderably different, with K.K. MC implementing a YFS-exponentiated algorithm, and PHOKHARA 
unexponentiated. Some comparisons of the MC's have been reported. 0, 0] Here, we will not 
compare the MC's, since that would test not just the comparison of the matrix elements, but also 
be affected by the presence or absence of exponentiation. Instead, we will implement both virtual 
photon corrections in the same MC, with the sole purpose of comparing the two matrix elements. 

Since the results of JMWY were developed for use in the YFS-based Monte Carlo (MC) program 
KK MC 0] ' we Wln look a t now the results of /C/C MC would change if the virtual correction of KR are 
substituted for those of JMWY in the same program. We calculate the YFS residuals [i\ as defined 

in refs. 0,0. The residuals of interest here are for single photon emission, ^ for i = 1,2 which 
were compared previously for several different virtual photon corrections in refs. 0, 0, 0, 0|. 
The superscript (i) denotes the power of a relative to the born level, (3^ = daBom/dQ. Thus, i = 1 
corresponds to tree-level hard photon emission, and i — 2 corresponds to the virtual corrections. The 
residuals are integrated up to a cutoff u ma x on the photon energy fraction to produce cross-sections 
denoted erv. 

Results for comparisons of these YFS residuals were previously reported in refs. @>EIE1>0>0] 
for the virtual photon corrections of refs. [l|,|3 at a CMS energy of 200 GeV. Some results were also 
reported in ref. |15| for 500 GeV. It was found in ref. that at 200 GeV, all of the fully differential 
results (excluding ref. 0, which is not fully differential) agreed to within 5 x 10~ 5 units of the Born 
cross section across the full range of v investigated (up to w max = 0.9999). Agreement to within 
2 x 10~ 4 units of the Born cross section was found at 500 GeV. 

Here, we will restrict our comparisons to the two fully differential results including mass cor- 
rections, JMWY and KR, discussed in the present paper. In addition to a comparison of the full 
results, we will compare the "massless" limits in which explicit mass terms in the virtual photon 
correction are omitted in both expressions. This will be useful in judging the size of the contribution 
of the mass corrections alone to the difference in the cross sections. Comparisons are made at CMS 
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FIG. 2: The differences between the exact results and the NLL result are shown at a CMS energy of 1 GeV 
both for (a) a cross section integrated to a cutoff u max and for (b) a partially differential cross section in v. 



energies of 1 GeV and 500 GeV. The 1 GeV scale was chosen as representative of some of the radia- 
tive return experiments for which PHOKHARA was designed, while the 500 GeV scale anticipates 
applications at a future linear collider (ILC). The 1 GeV comparisons are essentially new, since the 
only previously reported comparisons at low energy |18| were made before the matrix elements were 
completely implemented in their final form: it was found that some changes were needed to render 
the results sufficiently stable numerically when calculated for lower energy scales. The expressions 
shown in the previous sections are adequate for these calculations. 

Figure 1 shows the virtual correction cri^ — cri 1 ^, which is the pure 0(a 2 ) contribution to the 
single hard photon cross section. The results are obtained from 10 8 events created using the YFS3ff 
generator, which is the EEX3 option in the ICK, MC. The exact expressions are compared to a NLL 
approximation obtained from the expression eq. (|22[1 . All cross sections in these comparisons are 
normalized by dividing by the the radiationless Born cross-section OBorn for e + e~ — > It is 

seen that at both 1 GeV and 500 GeV, the agreement is very close between all the results compared. 

To see the differences between the various matrix elements, it is necessary to look at the NNLL 
contribution, since we know from the analytical results in the previous sections that these expressions 
agree to NLL order. Figure 2 examines the differences between the JMWY and KR results and the 
NLL contribution of eq. I|22|l at 1 GeV. This is a constant contribution (NNLL) to the cross section 
at 0(a 2 ) with no factors of the big logarithm L = \n(s/m 2 ). In FIG. 2(a), the difference between 
the integrated YFS residuals is shown, while in FIG. 2(b), the differences da^/dv are plotted. In 
other words, FIG. 2(b) is the derivative of the plot in FIG. 2(a), showing the difference between the 
virtual corrections at each photon energy v. It is found that the complete results of JMWY and KR 
differ by at most 2.2 x 10 -5 units of the Born cross section when integrated over the full range of v. 

The differential plot shows that the difference in the complete virtual correction is essentially 
zero at both endpoints v — 0, 1 and reaches a maximum of 3.3 x 10 -5 units of the Born cross section 
for v w 0.8. Most of the difference is found to be due to the mass correction terms. Without these 
terms, the results would agree to within at most 10 -5 , and much less over most of the range. The 
mass corrections agree in both the soft and hard limits at 1 GeV. The soft agreement mirrors that 
found analytically in the previous section. 

Figure 3 examines the NNLL contribution to the JMWY and KR results at 500 GeV. FIG. 3(a) 
shows the integrated distributions with the NLL contribution from eq. (|22|l subtracted. It is seen 
that the complete results, including mass corrections, differ by at most 1.6 x 10~ 5 units of the Born 
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FIG. 3: The differences between the exact results and the NLL result are shown at a CMS energy of 500 
GeV both for (a) a cross section integrated to a cutoff « ma x and for (b) a partially differential cross section 
in v. 



Cross section. Again, the bulk of the difference is seen to be due to the mass corrections. Without 
these, the results would agree to within 10~ 6 , except for the last bin. The last bin shows that the 
integrated contribution of the virtual photon factors over the entire range of photon energies differs 
by 5 x 1CP 6 units of the Born cross section. 

Differential results are compared in FIG. 3(b) up to v — 0.975. The final data point is omitted 
to permit the rest of the plot to be shown in better resolution. For the last data point (v — 0.9875), 
the KR result was 7.0 x 10~ 4 for the complete expression and 6.6 x 10" 4 for the "massless" part, 
while the JMWY result was 11.4 x 10~ 4 for the complete expression and 11.3 x 10 -4 for the massless 
part. The difference between the differential results at this point is likely to be enhanced due the 
steepness of the distribution. For the remaining points, the difference is less than 3 x 10~ 5 . In 
the limit when v — > 1, additional numerical issues may arise which are separate from the collinear 
singularities, so it is possible that some of the difference in this limit may be numerical. 

The absolute sizes of the NNLL contributions to the integals in FIG. 2(a) and FIG. 3(a) are 
also of interest, since these determine the range of v for which the much simpler NLL expression can 
be substituted for the complete result. For the JMWY results at 500 GeV, the NNLL contribution 
does not exceed 4 x 10~ 6 units of the Born cross section for cuts up to v max = 0.95, while for larger 
v, the NLL contribution can reach the 3.2 x 10 -5 level. At 1 GeV, the NNLL contribution to the 
JMWY result does not exceed 7 x 10~ 6 for the full range of photon energies. 

The fact that the NNLL contributions to the KR result are typically larger is mostly due to 
the mass corrections, which are considerably greater in the KR expression over most of the range of 
v. The mass corrections in the JMWY expression reach a level of at most 10~ 6 at 1 GeV, and are 
smaller at 500 GeV. In contrast, the KR mass corrections are typically of the order of 2 x 10~ 5 at 
intermediate photon energies. 

While we have expressed the comparisons in units of the non-radiative Born cross section, it is 
also useful to express them in units of the integrated cross section a\ SR for single real initial-state 
radiation, since there is always at least one photon in the radiative return context. Integrating the 
differential cross section of ref. 25], including the virtual photon contribution which cancels the 
IR divergence in the real emission cross section, gives a\ SR — 0.113 (JBorn at y/s — 1 GeV and 
crp R = 0.980 (JBorn at y/s = 500 GeV. Dividing by a^ R instead of <7Born would give essentially the 
same comparisons at 500 GeV, but increase the differences by a factor on the order of 10 at 1 GeV. 
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VI. CONCLUSIONS 



In this paper, we have compared two versions of the exact virtual corrections to initial-state 
bremsstrahlung in fermion pair production in detail, examining the collinear limits (NLL behavior) 
analytically both with and without the explicit mass corrections in each case. We found that the 
results of JMWY and KR agree to NLL order, both with and without masses, and that both results 
have the same soft limits. 

Numerical comparisons of the two virtual photon expressions were made by integrating them 
with the K.K. MC. The results were presented both in terms of integrated YFS residuals, to match 
comparisons made earlier, and in differential form, to permit comparison of the virtual photon 
corrections as a function of photon energy. It was found that the differential results agree to within 
3 x 10~ 5 units of the Born cross section across almost all photon energies, for both high energy (500 
GeV) and lower energy (1 GeV) scattering. The difference between the integrals of the YFS beta 
functions over the full photon energy range are also on the order of 10 -5 . This level of agreement is 
compatible with that found in earlier studies at 200 GeV. [THl flfi! UtI ITsj . 

The two virtual photon expressions compared agree much better when the mass corrections are 
omitted from both. This is because the mass corrections in the KR expression are typically of order 
10~ 5 , while the JMWY mass corrections tend to be of order 10 -6 units of the Born cross section 
over most of the photon energy range. Part of this difference could be due to the structure of the 
expressions, since the JMWY expression was created with the expectation of using it in the K.K. MC, 
while the KR result is foreign to this MC. However, some difference is not surprising, due to the 
difference in approach to the mass corrections. The JMWY result follows the approach of Berends 
et al [l^, adding the essential mass corrections for the collinear limits, and KR use an expansion in 
powers of m^/pi ■ k. The results are very different analytically, with the KR expressions containing 
many more terms. 

The results found in this paper give a clear estimate of the size of the difference in two exact, 
mass-corrected matrix elements for initial-state radiation in fermion pair production. We have also 
found that these matrix elements agree to NLL order, and have compared two approaches to adding 
mass corrections to the matrix element. Such information is important for estimating the precision 
of these matrix elements in applications to such processes as the LEP2 final data analyis, radiative 
return at «3> and B factories, and anticipated future ILC physics. 



APPENDIX A: STABILIZED FUNCTIONS 



The stabilized expressions in this paper depend on two sets of functions introduced to han- 
dle cancellations in differences of logarithms and dilogarithms with very similar arguments. The 
functions Lf n (x) and Sf n (x,y) are defined recursively by 

LfoOr) =ln(l + a;), U n+1 (x) = - (Lf n (x) - Lf n (0)) , (Ala) 

x 

Sf {x,y) = Sp(x + y), Sf n+1 (z,y) = - [Si n {x,y) - Sf^s.O)) ■ (Alb) 

V 

with Sp (x) the Spence dilogarithm function, which is also denoted Li2(x). 

Expansions can be used to evaluate Lf„(x) for x small or Si n (x,y) for y small. For the loga- 
rithmic difference function Lf„, 

°° (-~\\ k 

u n (x) = -J2 L jr-x k ~ n - (A2) 

k—n 

Differences between the functions Li n (x) for similar arguments may be calculated using the identity 
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fe=i 

valid for n > 1. Some identities for arguments containing ratios are 



- { Lf «(- + v) Lf «(-)i = YT^ Lfl {ttz) " 5> + y)fe " Lffc(a;) r (A3) 



+ = ,L fl (-^_), (A4a) 

(x + yf Lf 2 ( V -) = -x 2 Lf 2 f — O - x(x + y), (A4b) 
\ x / \ x I / 

(x + y) 3 Lf 3 (*) = x 3 Lf 3 (^-) + \xy{x + y). (A4c) 
\x / \x -\- y j z 

For the dilogarithmic difference function Sf„, 

00 k-n 

Sf n ( a; ,y) = ^LspW( a; ), (A5) 

where Sp ( fe '(a;) is the k th derivative of Sp (a;), which may be calculated recursively using 
Sp«(x) = Lfi(-aO, 
Sp (" +1 )(a;) = i | ["J J)n _ nS P (n) ( x )} for n ^ L ( A6 ) 
In the limit where both x and y are small, a double expansion is useful, 

Sf„(^ 2 /) = EE^(")- fc -V-"- (A7) 

The functions Lf„(x) have a logarithmic singularity at a; = —1. The functions Si n (x,y) are 
singular at x = 1 for n > 1, but can be calculated in for a; approaching 1 with x + y < 1 using the 
identity 

Sf n (x,y) = (-i)«+iSf n (l-a5,-y) + -^ T Lfi(-a;)Lf n (^ 



- J^T Lf.(x - 1) «„ (^) - Lf p (|) Lf„ ( 



(A8) 



where n > 1, p is the integer part of (n + l)/2 and g = 1 for n even, q = for n odd. A useful 
special case is 

Sf i (x,y) = S{ 1 (l-x, -y) - Lf 1 (x - 1) Lf j (-O - - ln(l - z - y) Lf i (¥■) . (A9) 
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